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Abstract 

Time-dependent wave-packet -propagation calculations are reported for the isotopic reactions, 
JJD- 1- CH 3 and D 2 + CH 3 , in six degrees of freedom and for zero total angular momentum. Initial 
state selected reaction probabilities for different initial rotational-vibrational states are presented 
in this study. This study shows that excitations of the HD{D 2 ) enhances the reactivities; whereas 
the excitations of the CH 3 umbrella mode have the opposite effects. This is consistent with the re- 
action of H 2 + CH 3 . The comparison of these three isotopic reactions also shows the isotopic effects 
in the initial-state-selected reaction probabilites. The cumulative reaction probabilities (CRP) are 
obtained by summing over initial-state-selected reaction probabilities. The energy-shift approx- 
imation to account for the contribution of degrees of freedom missing in the six dimensionality 
calculation is employed to obta n approximate full-dimensional CRPs. The rate constant compar- 
ison shows H 2 + CH 3 reaction has the biggest reactivity, then HD + CH 3 , and D 2 + CH 3 has the 
smallest. 
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I. INTRODUCTION 


Due to the importance of the reaction in methane combustion and pyrolysisfl], H + 
CH a -4 H 2 + CH 3 , and its reverse reaction, H 2 + CH 3 -4 i/ + CH 4 have been studied 
experimentally[2-l 1] and theoretically[13-18] since the 1960s. In addition, the reverse reac- 
tion also plays an important role in the modeling of the chemistry of hydrocarbons in the 
atmospheres of giant planets[12]. Besides a number of quantum dynamics calculations [19— 
26] in reduced dimensionality and a full-dimensional transition state flux method quantum 
calculation on the forward leaction H + CH4, recently we reported a 6D time-dependent 
quantum study of the reverse reaction H 2 + CH 3 -4 H + CH 4 [ 27], 

In this reduced dimensional calculation, CH 3 was regarded as one pseudo-diatom 
molecule; wave-packet propagation method has been employed to calculate the reaction 
probabilities; energy shifting[28-30] methods are used to obtain approximate full dimen- 
sionality cumulative reaction probability for zero angular momentum J, and J-K shifting 
methods[28-30] are employed to calculate total cumulative reaction probability for the full 
dimensionality rate constants. Initial state selected reaction probabilities were calculated. 
We found that the excitation of the H 2 substantially enhances the reactivity, whereas the 
excitation of the umbrella n otion has the opposite effects. We also found that the ji=0 
initial rotational state of H 2 gives the largest reaction probability of the rotational excita- 
tion states of H 2 . On the otaer hand, j 2 =l initial rotational state of CH 3 has the largest 
reaction probability of the CH 3 rotational states. Rate constant from our reduced dimen- 
sional calculation is larger than currently available experimental data, especially at the lower 
temperature range. 

The quantum calculation mentioned above was done using the twelve degrees of freedom 
potential energy surface of Jordan and Gilbert [31]. This potential has a C 3v transition state 
with a collinear structure of H — H — C and the center of mass of H 3 . The vibrationally 
adiabatic barrier height of tl is potential energy surface is 10.9 kcal/mol for the forward 
reaction, H + CH 4 -4 H 2 + CH 3 and 11.0 kcal/mol for the reverse reaction, H 2 + CH 3 — > 
H + CH 4 , respectively. This potential energy surface is based in part on previous ab initio 
calculations of the saddle point[l7], as well as on a previous, semiglobal potential used in 
variational transition state theory calculations of the rate constant[32]. The rate constant 
from the theoretical calculations of both the forward react,ion[19, 20, 22-25] and the reverse 
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reaction[27], comparing to the experiments, indicates that the barrier height on the Jordan- 
Gilbert potential may be to ) low. 

In this paper, we extend the previous quantum dynamics study of the reaction H 2 +CH 3 — > 
H + CHi [27]to include the isotopic reactions HD + CH$ and D 2 + CH S , by treating them 
as diatom-diatom reaction systems and carry out full dimensional quantum dynamics calcu- 
lations for those pseudo diatom-diatom reactions. Initial state selected reaction probability, 
cumulative reaction probability(CRP), and thermal rate constant are calculated. Ground 
state reaction probabilities, GRP, rate constants are compared in this report for these three 
reactions. 

This paper is organized is follows: Section II briefly reviews the theory used in this 
calculation and the energy-shifting methods are also presented. Section III gives numerical 
results of the calculation, and section IV concludes. 

II. THEORY AND METHODS 

A. Time-dependent wive-packet-propagation method 

This study employs the time-dependent wave-packet propagation method[27, 33-35] to 
carry out the quantum dynamics calculation. The two reactions, HD + CH$ and D 2 + CH S 
are treated as diatom-diatom reactions HD(D 2 ) -I- CX, as it is defined in the study of the 
reaction H 2 + CH 3 [27], where X is a pseudodiatom with the mass of the three nonreactive 
hydrogen atoms and located at their center of mass. The vibrational motion of CX approx- 
imately represents the umbrella motion of CH 3 , where the three hydrogen atoms move in 
phase , and the C — H bone length in CH 3 is fixed at its equilibrium geometries. Thus 
the reduced dimensionality cilculation for this pseudo diatom-diatom uses six degrees of 
freedom. 

The mathematical treatment for these tw r o reactions is the same as for the previous 
calculations of H 2 + CH 3 . A direct product of a specific initial rotational-vibrational eigen- 
functions of HD(D 2 ) and CX and a localized standard Gaussian function is chosen as the 
initial wave-packet. We employ the split-operator method to propagate the time-dependent 
-wavefuntion. And the initial-s tate-selected reaction probability is extracted from the time- 
dependent reactive flux in the standard way. 
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B. Cumulative reaction probability and thermal rate constant 

In order to obtain an approximate full dimensional (12D) cumulative reaction probabil- 
ity(CRP) for J=0 and thermal rate constant, we first calculate the reduced dimensional(6D) 
CRP(J=0), denoted N J=0 ’ rc (E), by summing over all the initial-state-selected rovibrational 
reaction probability P^ nr (E), 

N' =0,Td (E) = Y, £ p ^°n,SE) (1) 

"1 ,l/ 2 =0 ji,h=0 

where E is the total Energy. 

Our reduced dynamics calculation has six degrees of freedom, while the title reaction 
system has twelve degrees o: freedom. So the energy shifting method has to be applied to 
obtain the CRP for the full c dimensionality using our reduced dimensional calculation. Thus 


N J=o,m, = V „). (2) 

here v \V 2 ,...,vq are the quantum numbers of the six vibrational degrees of freedom of the 
transition state, which are not included in our reduced dimensional calculation. These 
energies are given in the harmonic approximation by 

i/e = Y fa + b- (3) 

t = l 

The J-K shifting method is employed here to calculate the CRP for nonzero total angular 
momentum J. The total CRP. N(E) is defined as the sum of CRP of all the open J and K 
channels 


N(E) = £(2 J + 1) Y N J ^“'(B - E< jk ), (4) 

7-0 K=-J 

where E^ JK is the rigid rotor rc tational energy of the reaction system at the transition state. 
In the present case, the system at transition state is a well defined symmetric top 
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£jk — B*J(J + 1) + {A x - B X )K 2 . (5) 

A 1 and B x are the rotationa constants of H 2 CH 3 at the full dimensional transition state. 

The thermal rate constant can be computed exactly from the expression 

1 r°° 

k(T) = — / N(E)e- E ' kBT dE. (6) 

react Jo 

The above equation simplifies in the J-K shifting approximation, Eq. (4), to 

O x r°° 

k(T) =■■ - ~ - TOt / N(E) J=0JuU e- E/kBT dE, (7) 

react J 0 

where Q x rot is the rotational partition function of the reaction system at the transition state 
and Q react is the reactant partition function, which is written as a product of vibrational, 
rotational, and translational partition functions. As noted in Ref. [24], the consistent 
approach to calculate the vibrational partition function of CH 3 is to combine the results 
of the reduced dimensional quantum calculations with the ’missing’ vibrational frequencies 
which are not included in the reduced dimensionality. The purpose of doing this is to be 
consistent with the procedure that we use to convert the full dimensional CRP from the 
reduced dimensional CRP in Eq. (2). 

III. RESULTS AND DISCUSSION 
A. Numerical aspects 

Sixty Sine basis functions are used to expand the translational coordinates, 30 potential 
optimized vibrational functions for the HD(D 2 ) vibrational coordinates, 8 potential opti- 
mized vibrational functions fir the CX vibrational coordinates and and total 1092 coupled 
angular momentum basis usee for the angular part. The rotational constants, A x and B x 
in Eq. (5), are 5.12 cm -1 and 1.33 cm -1 for HD + CH 3 system; 5.12 cm -1 and 1.23 cm~ l 
for D 2 T CH 3 system. The \ibrational frequencies that are ’missing’ in our reduced di- 
mensional calculations at the transition state are: 589 cm -1 (2), 1439 cm' 1 ] 2), 2959cm -1 
and 3092cm -1 (2)[24]. These missing’ vibrational frequencies are used to obtain the full 
dimensionality CRP for J=0 ir Eq. (2). 
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B. Initial-state-selected probability 


Fig.l shows the reaction probabilities for the first four initial vibrational states of H D{u i, 
ji=0) with CX at ground state as a function of kinetic energy. As seen, vibrational excitation 
of HD enhances this reactivity with 1^=3 having the largest reaction probability. And in 
addition, vibrational excitation of HD also enhances the reactivity by lowing the reaction 
threshold. The same trend of reaction probabilities is also shown in Fig. 2 for reaction 
;'i=0) + CX (j / 2=0, 

The reaction probabilities of vibrational excitation of CX(u 2 , .72=0, ^=^=0) are plotted 
in Fig.3 and Fig.4 respectively for HD + CH 3 and D 2 + CH 3 . Comparing to Fig.l, These 
figures simply tell us that excitations of the umbrella motion of CX hinder the reactivity 
due to CX mainly functions as a spector in this reaction. The more of the CH 3 vibrates, 
the harder for atoms in HD and D 2 to ’catch’ the C atom in CH 3 . 

The comparison of reaction probabilities for the three reactions, H 2 (KD,D 2 ) + CH$ with 
all reactants at their ground ro-vibrational states are plotted in Fig. 5. The isotopic effect 
has been shown in this figure that the H 2 + CH 3 reaction has the smallest threshold, and 
D 2 + CH 3 has the highest. One can say, the H 2 + C H 3 reaction has the lowest adiabatic 
barrier because of the largest zero-point energy and D 2 + CH 3 reaction has higest adiabatic 
barrier. So at relatively lower kinetic energies, H 2 + CH Z reaction has the largest reaction 
probability, D 2 + CH 3 reaction has the smallest reaction probability. However, the reaction 
probablities have the reverse order at relatively high kinetic energies. This is also caused by 
the isotopic effect because D - is the heaviest among the three hydorgen isotopes. And at 
relatively high kinetic energies, it moves slower than the other two, H 2 and HD, respect to 
CH 3 , so it’s easiler to get combined with the C atom in CH 3 to let reaction occur. 

C. Cumulative reaction probability and thermal rate constant 

A normal mode analysis is performed to determine the ’missing’ modes in our reduced 
dimensional quantum calculation, as also seen in Ref. [24], Then Eq.(2) is used to obtain 
the full dimensional CRP for , =0. For total scattering energy up to 2.0 eV, all the missing 
modes given in Section A. contribute to convert the full dimensional CRP from the six- 
degree-of-freedom quantum calculation. And the comparison of full dimensional converted 
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CPR(J=0) is shown in Fig. 6 as a function of the total energy, which includes the energy 
difference 0.702 eV which not included in the 6D reduced dimensional ZPE. Using the 
harmonic approximation of ihe ’missing’ degrees of freedom in CH 3 is consistent with using 
this approximation in the energy-shifting expression in Eq. (2). This figure shows that 
D 2 + CH$ has the largest C RP(J=0), which is expected since D 2 + CH 3 reaction has more 
open channels (initial rotational-vibrational states) than the other two in the energy rangy 
we are interested here. 

Our previous study[27] showed the rate constant of H 2 + CH 3 is larger the experiment 
rates, especially in the low temperature, this indicates the barrier on the Jordan-Gilbert 
potential is too low. Compaiision of our present rate constants of the three reactions shows 
H 2 + CH 3 the largest rate constant, then HD + CH 3 and Do + CH 3 the smallest in Fig. 7. 
In other words, the isotopic effect shows, from these three reactions, that the larger of the 
mass of the reactant particle the smaller reactivity. 

IV. SUMMARY AND CONCLUSIONS 

In this study, we perform the isotopic chemical reaction study of HD + CH 2 and D 2 + 
CH 3 , which is an extension study of the H 2 + CH 2 reaction. Time-dependent wave-packet 
propagation method is carried out in six-degree-of freedom for the title reactions. The 
energy-shifting method has been employed to obtain the full dimensional cumulative reaction 
probability for J=0. And J-K shifting approach is applied to this CRP to obtain the thermal 
rate constant which is compared with different experiments. A potential surface due to 
Jordan and Gilbert, which h;is been used in previous quantum calculation of the forward 
reaction H + CH 4 — r H 2 + C H 3 , -was used in these calculations. 

Initial-state-selected reaction probability and the cumulative reaction probability were 
calculated for zero total angular momentum. The vibrational excitations of HD(D 2 ) sub- 
stantially enhance the reactio i probabilities. On the contrary, the vibrational excitation of 
CH 3 causes the opposite effects in these reactions. The H 2 + CH 2 reaction has the smallest 
adiabatic barrier due to the largest ZPE because of the isotopic effect. This study also shows 
at relatively high kinetic energy, D 2 + CH :i is the most reactive one. 

An energy shifting method was used to convert the reduced dimensional CRP for zero 
total angular momentum to a full dimensional one. Our rate constant comparison shows 
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H 2 + CH$ reaction has the largest rate constant against the smallest of D 2 + CH 3 . 
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Fig. 1 The initial state selected reaction probability for the reaction HD(i/\,ji = 0) + 
C H 3 {v 2 =j2=0) as a function of kinetic energy (J=0). 

Fig. 2 Same as Fig. 1 but fo: the reaction ^ 2(^1 Ji=0)+CH 3 (u2=j2=0). 

Fig. 3 The initial state selected reaction probability for the reaction FfD(n' 1 =0ji=0) + 
CH 3 (v u j 2 =0) as a furction of kinetic energy (J=0). 

Fig. 4 Same as Fig. 3 but foi the reaction CH z {y x j 2 =0). 

Fig. 5 Comparison of the ground state reaction probability for all three isotopic reaction: 
H 2 {HD , Z? 2 ) + CH Z reaction as a function of kinetic energy. 

Fig. 6 Comparison of cumulative reaction probabilities N J=0 ^ ull (E) for all three isotopic 
reactions as a function of total energy. 

Fig. 7 Comparison of all three isotopic thermal rate constant. 
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